Parallel evolution despite low genetic diversity in three-spined sticklebacks.

When populations repeatedly adapt to similar environments they can evolve similar phenotypes based on shared genetic mechanisms (parallel evolution). The likelihood of parallel evolution is affected by demographic history, as it depends on the standing genetic variation of the source population. The three-spined stickleback (Gasterosteus aculeatus) repeatedly colonized and adapted to brackish and freshwater. Most parallel evolution studies in G. aculeatus were conducted at high latitudes, where freshwater populations maintain connectivity to the source marine populations. Here, we analysed southern and northern European marine and freshwater populations to test two hypotheses. First, that southern European freshwater populations (which currently lack connection to marine populations) lost genetic diversity due to bottlenecks and inbreeding compared to their northern counterparts. Second, that the degree of genetic parallelism is higher among northern than southern European freshwater populations, as the latter have been subjected to strong drift due to isolation. The results show that southern populations exhibit lower genetic diversity but a higher degree of genetic parallelism than northern populations. Hence, they confirm the hypothesis that southern populations have lost genetic diversity, but this loss probably happened after they had already adapted to freshwater conditions, explaining the high degree of genetic parallelism in the south.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord.

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

Condition- and context-dependent variation of sexual dimorphism across lizard populations at different spatial scales.

The evolution of sexual dimorphism (SD) is driven by intricate interplays between sexual and natural selection. When it comes to SD variation within populations, however, environmental factors play a major role. Sexually selected traits are expected to be strongly dependent on individual body condition, which is influenced by the local environment that individuals experience. As a consequence, the degree of SD may also depend on resource availability. Here, we investigated the potential drivers of SD expression at two sexually dimorphic morphometric traits, body size (snout vent length) and head shape (head geometric morphometrics), in the Italian wall lizard (Podarcis siculus). We assessed the existence of condition- and context-dependent SD across ten islands of the Aeolian archipelago (southern Italy), at within- and among-population scales. We observed strong geographical variation of SD among islands, and tested three potential SD predictors related to resource availability (individual body condition, ecosystem productivity, temperature). Body condition and ecosystem productivity were the main drivers of body size SD variation, and body condition was also the main driver for head shape SD. Our results highlight that the expression of SD in the Italian wall lizard is both condition- and context-dependent. These results are congruent at within- and among-populations scales highlighting that spatial multi-scale analysis represents a useful approach to understand patterns of SD expression.

Allopatric origin of sympatric whitefish morphs with insights on the genetic basis of their reproductive isolation.

The European whitefish (Coregonus lavaretus) species complex is a classic example of recent adaptive radiation. Here, we examine a whitefish population introduced to northern Finnish Lake Tsahkal in the late 1960s, where three divergent morphs (viz. littoral, pelagic, and profundal feeders) were found 10 generations after. Using demographic modeling based on genomic data, we show that whitefish morphs evolved during a phase of strict isolation, refuting a rapid sympatric divergence scenario. The lake is now an artificial hybrid zone between morphs originated in allopatry. Despite their current syntopy, clear genetic differentiation remains between two of the three morphs. Using admixture mapping, we identify five SNPs associated with gonad weight variation, a proxy for sexual maturity and spawning time. We suggest that ecological adaptations in spawning time evolved in allopatry are currently maintaining partial reproductive isolation in the absence of other barriers to gene flow.

Genomic insights into the historical and contemporary demographics of the grey reef shark.

Analyses of genetic diversity can shed light on both the origins of biodiversity hotspots, as well as the conservation status of species that are impacted by human activities. With these objectives, we assembled a genomic dataset of 14,935 single nucleotide polymorphisms from 513 grey reef sharks (Carcharhinus amblyrhynchos) sampled across 17 locations in the tropical Indo-Pacific. We analysed geographic variation in genetic diversity, estimated ancient and contemporary effective population size (Ne) across sampling locations (using coalescent and linkage disequilibrium methods) and modelled the history of gene flow between the Coral Triangle and the Coral Sea. Genetic diversity decreased with distance away from the Coral Triangle and north-western Australia, implying that C. amblyrhynchos may have originated in this region. Increases in Ne were detected across almost all sampling locations 40,000-90,000 generations ago (approximately 0.6-1.5 mya, given an estimated generation time of 16.4 years), suggesting a range expansion around this time. More recent, secondary increases in Ne were inferred for the Misool and North Great Barrier Reef sampling locations, but joint modelling did not clarify whether these were due to population growth, migration, or both. Despite the greater genetic diversity and ancient Ne observed at sites around Australia and the Coral Triangle, remote reefs around north-western New Caledonia had the highest contemporary Ne, demonstrating the importance of using multiple population size assessment methods. This study provides insight into both the past and present demographics of C. amblyrhynchos and contributes to our understanding of evolution in marine biodiversity hotspots.

Cast away in the Adriatic: Low degree of parallel genetic differentiation in three-spined sticklebacks.

The three-spined stickleback (Gasterosteus aculeatus) has repeatedly and independently adapted to freshwater habitats from standing genetic variation (SGV) following colonization from the sea. However, in the Mediterranean Sea G. aculeatus is believed to have gone extinct, and thus the spread of locally adapted alleles between different freshwater populations via the sea since then has been highly unlikely. This is expected to limit parallel evolution, that is the extent to which phylogenetically related alleles can be shared among independently colonized freshwater populations. Using whole genome and 2b-RAD sequencing data, we compared levels of genetic differentiation and genetic parallelism of 15 Adriatic stickleback populations to 19 Pacific, Atlantic and Caspian populations, where gene flow between freshwater populations across extant marine populations is still possible. Our findings support previous studies suggesting that Adriatic populations are highly differentiated (average FST  ≈ 0.45), of low genetic diversity and connectivity, and likely to stem from multiple independent colonizations during the Pleistocene. Linkage disequilibrium network analyses in combination with linear mixed models nevertheless revealed several parallel marine-freshwater differentiated genomic regions, although still not to the extent observed elsewhere in the world. We hypothesize that current levels of genetic parallelism in the Adriatic lineages are a relic of freshwater adaptation from SGV prior to the extinction of marine sticklebacks in the Mediterranean that has persisted despite substantial genetic drift experienced by the Adriatic stickleback isolates.

Population Structure Limits Parallel Evolution in Sticklebacks.

Population genetic theory predicts that small effective population sizes (Ne) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e., parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically codistributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbors less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius: 2,996 versus 379 single nucleotide polymorphisms located within 26 versus 9 genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius, respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution.

Biases in Demographic Modeling Affect Our Understanding of Recent Divergence.

Testing among competing demographic models of divergence has become an important component of evolutionary research in model and non-model organisms. However, the effect of unaccounted demographic events on model choice and parameter estimation remains largely unexplored. Using extensive simulations, we demonstrate that under realistic divergence scenarios, failure to account for population size (Ne) changes in daughter and ancestral populations leads to strong biases in divergence time estimates as well as model choice. We illustrate these issues reconstructing the recent demographic history of North Sea and Baltic Sea turbots (Scophthalmus maximus) by testing 16 isolation with migration (IM) and 16 secondary contact (SC) scenarios, modeling changes in Ne as well as the effects of linked selection and barrier loci. Failure to account for changes in Ne resulted in selecting SC models with long periods of strict isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for Ne changes suggest recent (<6 kya) divergence with constant gene flow. We further show how interpreting genomic landscapes of differentiation can help discerning among competing models. For example, in the turbot data, islands of differentiation show signatures of recent selective sweeps, rather than old divergence resisting secondary introgression. The results have broad implications for the study of population divergence by highlighting the potential effects of unmodeled changes in Ne on demographic inference. Tested models should aim at representing realistic divergence scenarios for the target taxa, and extreme caution should always be exercised when interpreting results of demographic modeling.

The roles of climate, geography and natural selection as drivers of genetic and phenotypic differentiation in a widespread amphibian Hyla annectans (Anura: Hylidae).

The role of geological events and Pleistocene climatic fluctuations as drivers of current patterns of genetic variation in extant species has been a topic of continued interest among evolutionary biologists. Nevertheless, comprehensive studies of widely distributed species are still rare, especially from Asia. Using geographically extensive sampling of many individuals and a large number of nuclear single nucleotide polymorphisms (SNPs), we studied the phylogeography and historical demography of Hyla annectans populations in southern China. Thirty-five sampled populations were grouped into seven clearly defined genetic clusters that closely match phenotype-based subspecies classification. These lineages diverged 2.32-5.23 million years ago (Ma), a timing that closely aligns with the rapid and drastic uplifting of the Qinghai-Tibet Plateau and adjacent southwest China. Demographic analyses and species distribution models indicate that different populations of this species have responded differently to past climatic changes. In the Hengduan Mountains, most populations experienced a bottleneck, whereas the populations located outside of the Hengduan Mountains have gradually declined in size since the end of the last glaciation. In addition, the levels of phenotypic and genetic divergence were strongly correlated across major clades. These results highlight the combined effects of geological events and past climatic fluctuations, as well as natural selection, as drivers of contemporary patterns of genetic and phenotypic variation in a widely distributed anuran in Asia.

On the causes of geographically heterogeneous parallel evolution in sticklebacks.

The three-spined stickleback (Gasterosteus aculeatus) is an important model system for the study of parallel evolution in the wild, having repeatedly colonized and adapted to freshwater from the sea throughout the northern hemisphere. Previous studies identified numerous genomic regions showing consistent genetic differentiation between freshwater and marine ecotypes but these had typically limited geographic sampling and mostly focused on the Eastern Pacific region. We analysed population genomic data from global samples of the three-spined stickleback marine and freshwater ecotypes to detect loci involved in parallel evolution at different geographic scales. Most signatures of parallel evolution were unique to the Eastern Pacific and trans-oceanic marine-freshwater differentiation was restricted to a limited number of shared genomic regions, including three chromosomal inversions. On the basis of simulations and empirical data, we demonstrate that this could result from the stochastic loss of freshwater-adapted alleles during the invasion of the Atlantic basin and selection against freshwater-adapted variants in the sea, both of which can reduce standing genetic variation available for freshwater adaptation outside the Eastern Pacific region. Moreover, the elevated linkage disequilibrium associated with marine-freshwater differentiation in the Eastern Pacific is consistent with secondary contact between marine and freshwater populations that evolved in isolation from each other during past glacial periods. Thus, contrary to what earlier studies from the Eastern Pacific region have led us to believe, parallel marine-freshwater differentiation in sticklebacks is far less prevalent and pronounced in all other parts of the species global distribution range.

The evolution of polymorphism in the warning coloration of the Amazonian poison frog Adelphobates galactonotus.

While intraspecific variation in aposematic signals can be selected for by different predatory responses, their evolution is also contingent on other processes shaping genetic variation. We evaluate the relative contributions of selection, geographic isolation, and random genetic drift to the evolution of aposematic color polymorphism in the poison frog Adelphobates galactonotus, distributed throughout eastern Brazilian Amazonia. Dorsal coloration was measured for 111 individuals and genetic data were obtained from 220 individuals at two mitochondrial genes (mtDNA) and 7963 Single Nucleotide Polymorphisms (SNPs). Four color categories were described (brown, blue, yellow, orange) and our models of frog and bird visual systems indicated that each color was distinguishable for these taxa. Using outlier and correlative analyses we found no compelling genetic evidence for color being under divergent selection. A time-calibrated mtDNA tree suggests that the present distribution of dorsal coloration resulted from processes occurring during the Pleistocene. Separate phylogenies based on SNPs and mtDNA resolved the same well supported clades, each containing different colored populations. Ancestral character state analysis provided some evidence for evolutionary transitions in color type. Genetic structure was more strongly associated with geographic features, than color category, suggesting that the distribution of color is explained by localized processes. Evidence for geographic isolation together with estimates of low effective population size implicates drift as playing a key role in color diversification. Our results highlight the relevance of considering the neutral processes involved with the evolution of traits with important fitness consequences.

Estimating uncertainty in divergence times among three-spined stickleback clades using the multispecies coalescent.

Incomplete lineage sorting (ILS) can lead to biased divergence time estimates. To explore if and how ILS has influenced the results of a recent worldwide phylogeny of three-spined sticklebacks (Gasterosteus aculeatus), we estimated divergence times among major clades by applying both a concatenation approach and the multispecies coalescent (MSC) model to single-nucleotide polymorphisms. To further test the influence of different calibration strategies, we applied different calibrations to the root and to younger nodes in addition to the ones used in the previous study. Both the updated calibrations and the application of the MSC model influenced divergence time estimates, sometimes significantly. The new divergence time estimates were more ancient than in the previous study for older nodes, whereas the estimates of younger nodes were not strongly affected by the re-analyses. However, since the applied MSC method employs a simple substitution model and cannot account for changes in population size, we suggest that different analytical approaches and calibration strategies should be used in order to explore uncertainty in divergence time estimates. This study provides a valuable reference timeline for the ages of worldwide three-spined stickleback populations and emphasizes the need to embrace, rather than obscure, uncertainties around divergence time estimates.

A phylogenomic perspective on diversity, hybridization and evolutionary affinities in the stickleback genus Pungitius.

Hybridization and convergent evolution are phenomena of broad interest in evolutionary biology, but their occurrence poses challenges for reconstructing evolutionary affinities among affected taxa. Sticklebacks in the genus Pungitius are a case in point: evolutionary relationships and taxonomic validity of different species and populations in this circumpolarly distributed species complex remain contentious due to convergent evolution of traits regarded as diagnostic in their taxonomy, and possibly also due to frequent hybridization among taxa. To clarify the evolutionary relationships among different Pungitius species and populations globally, as well as to study the prevalence and extent of introgression among recognized species, genomic data sets of both reference genome-anchored single nucleotide polymorphisms and de novo assembled RAD-tag loci were constructed with RAD-seq data. Both data sets yielded topologically identical and well-supported species trees. Incongruence between nuclear and mitochondrial DNA-based trees was found and suggested possibly frequent hybridization and mitogenome capture during the evolution of Pungitius sticklebacks. Further analyses revealed evidence for frequent nuclear genetic introgression among Pungitius species, although the estimated proportions of autosomal introgression were low. Apart from providing evidence for frequent hybridization, the results challenge earlier mitochondrial and morphology-based hypotheses regarding the number of species and their affinities in this genus: at least seven extant species can be recognized on the basis of genetic data. The results also shed new light on the biogeographical history of the Pungitius-complex, including suggestion of several trans-Arctic invasions of Europe from the Northern Pacific. The well-resolved phylogeny should facilitate the utility of this genus as a model system for future comparative evolutionary studies.

Cryptic temporal changes in stock composition explain the decline of a flounder (Platichthys spp.) assemblage.

Unobserved diversity, such as undetected genetic structure or the presence of cryptic species, is of concern for the conservation and management of global biodiversity in the face of threatening anthropogenic processes. For instance, unobserved diversity can lead to overestimation of maximum sustainable yields and therefore to overharvesting of the more vulnerable stock components within unrecognized mixed-stock fisheries. We used DNA from archival (otolith) samples to reconstruct the temporal (1976-2011) genetic makeup of two mixed-stock flounder fisheries in the Åland Sea (AS) and the Gulf of Finland (GoF). Both fisheries have hitherto been managed as a single stock of European flounders (Platichthys flesus), but were recently revealed to target two closely related species: the pelagic-spawning P. flesus and the newly described, demersal-spawning P. solemdali. While the AS and GoF fisheries were assumed to consist exclusively of P. solemdali, P. flesus dominated the GoF flounder assemblage (87% of total) in 1983, had disappeared (0%) by 1993, and remained in low proportions (10%-11%) thereafter. In the AS, P. solemdali dominated throughout the sampling period (>70%), and P. flesus remained in very low proportions after 1983. The disappearance of P. flesus from the GoF coincides in time with a dramatic (~60%) decline in commercial landings and worsening environmental conditions in P. flesus' northernmost spawning ground, the Eastern Gotland Basin, in the preceding 4-6 years. These results are compatible with the hypothesis that P. flesus in the GoF is a sink population relying on larval subsidies from southern spawning grounds and the cause of their disappearance is a cessation of larval supply. Our results highlight the importance of uncovering unobserved genetic diversity and studying spatiotemporal changes in the relative contribution of different stock components, as well as the underlying environmental causes, to manage marine resources in the age of rapid anthropogenic change.

Worldwide phylogeny of three-spined sticklebacks.

Stickleback fishes in the family Gasterosteidae have become model organisms in ecology and evolutionary biology. However, even in the case of the most widely studied species in this family - the three-spined stickleback (Gasterosteus aculeatus) - the worldwide phylogenetic relationships and colonization history of the different populations and lineages remain poorly resolved. Using a large collection of samples covering most parts of the species distribution range, we subjected thousands of single nucleotide polymorphisms to coalescent analyses in order to reconstruct a robust worldwide phylogeny of extant G. aculeatus populations, as well as their ancestral geographic distributions using Statistical-Dispersal Vicariance and Bayesian Binary MCMC analyses. The results suggest that contemporary populations originated from the Pacific Ocean in the Late Pleistocene, and the Atlantic was colonized through the Arctic Ocean by a lineage that diverged from Pacific sticklebacks ca 44.6 Kya. This lineage contains two branches: one that is distributed in the Mediterranean area, from the Iberian Peninsula to the Black Sea ('Southern European Clade'), and another that is comprised of populations from northern Europe and the east coast of North America ('Trans-Atlantic Clade'). Hence, the results suggest that the North American East Coast was colonized by trans-Atlantic migration. Coalescence-based divergence time estimates suggest that divergence among major clades is much more recent than previously estimated.

Extraordinarily rapid speciation in a marine fish.

Divergent selection may initiate ecological speciation extremely rapidly. How often and at what pace ecological speciation proceeds to yield strong reproductive isolation is more uncertain. Here, we document a case of extraordinarily rapid speciation associated with ecological selection in the postglacial Baltic Sea. European flounders (Platichthys flesus) in the Baltic exhibit two contrasting reproductive behaviors: pelagic and demersal spawning. Demersal spawning enables flounders to thrive in the low salinity of the Northern Baltic, where eggs cannot achieve neutral buoyancy. We show that demersal and pelagic flounders are a species pair arising from a recent event of speciation. Despite having a parapatric distribution with extensive overlap, the two species are reciprocally monophyletic and show strongly bimodal genotypic clustering and no evidence of contemporary migration, suggesting strong reproductive isolation. Divergence across the genome is weak but shows strong signatures of selection, a pattern suggestive of a recent ecological speciation event. We propose that spawning behavior in Baltic flounders is the trait under ecologically based selection causing reproductive isolation, directly implicating a process of ecological speciation. We evaluated different possible evolutionary scenarios under the approximate Bayesian computation framework and estimate that the speciation process started in allopatry ∼2,400 generations ago, following the colonization of the Baltic by the demersal lineage. This is faster than most known cases of ecological speciation and represents the most rapid event of speciation ever reported for any marine vertebrate.

Effects of temperature, salinity and composition of the dinoflagellate assemblage on the growth of Gambierdiscus carpenteri isolated from the Great Barrier Reef.

Increases in reported incidence of ciguatera fish poisoning (hereafter ciguatera) have been linked to warmer sea temperatures that are known to trigger coral bleaching events. The drivers that trigger blooms of ciguatera-causing dinoflagellates on the Great Barrier Reef (GBR) are poorly understood. This study investigated the effects of increased temperatures and lowered salinities, often associated with environmental disturbance events, on the population growth of two strains of the potentially ciguatera-causing dinoflagellate, Gambierdiscus carpenteri (NQAIF116 and NQAIF380). Both strains were isolated from the central GBR with NQAIF116 being an inshore strain and NQAIF380 an isolate from a stable environment of a large coral reef aquarium exhibit in ReefHQ, Townsville, Australia. Species of Gambierdiscus are often found as part of a mixed assemblage of benthic toxic dinoflagellates on macroalgal substrates. The effect of assemblage structure of dinoflagellates on the growth of Gambierdiscus populations has, however, not been explored. The study, therefore investigated the growth of G. carpenteri within mixed assemblages of benthic dinoflagellates. Population growth was monitored over a period of 28days under three salinities (16, 26 and 36) and three temperature (24, 28 and 34°C) conditions in a fully crossed experimental design. Temperature and salinity had a significant effect on population growth. Strain NQAIF380 exhibited significantly higher growth at 28°C compared to strain NQAIF116, which had highest growth at 24°C. When strain NQAIF116 was co-cultured with the benthic dinoflagellates, Prorocentrum lima and Ostreopsis sp., inhibitory effects on population growth were observed at a salinity of 36. In contrast, growth stimulation of G. carpenteri (strain NQAIF116) was observed at a salinity of 26 and particularly at 16 when co-cultured with Ostreopsis-dominated assemblages. Range expansion of ciguatera-causing dinoflagellates could lead to higher frequency of reported ciguatera illness in populated temperate Australian regions, outside the tropical range of the GBR. Therefore, the findings on salinity and temperature tolerance of two strains of G. carpenteri indicates potential adaptability to different local environmental conditions. These are baseline data for future investigations into the potential southward range expansion of ciguatera-causing dinoflagellates originating from the GBR.

Connectivity in grey reef sharks (Carcharhinus amblyrhynchos) determined using empirical and simulated genetic data.

Grey reef sharks (Carcharhinus amblyrhynchos) can be one of the numerically dominant high order predators on pristine coral reefs, yet their numbers have declined even in the highly regulated Australian Great Barrier Reef (GBR) Marine Park. Knowledge of both large scale and fine scale genetic connectivity of grey reef sharks is essential for their effective management, but no genetic data are yet available. We investigated grey reef shark genetic structure in the GBR across a 1200 km latitudinal gradient, comparing empirical data with models simulating different levels of migration. The empirical data did not reveal any genetic structuring along the entire latitudinal gradient sampled, suggesting regular widespread dispersal and gene flow of the species throughout most of the GBR. Our simulated datasets indicate that even with substantial migrations (up to 25% of individuals migrating between neighboring reefs) both large scale genetic structure and genotypic spatial autocorrelation at the reef scale were maintained. We suggest that present migration rates therefore exceed this level. These findings have important implications regarding the effectiveness of networks of spatially discontinuous Marine Protected Areas to protect reef sharks.

The diversity of Coolia spp. (Dinophyceae Ostreopsidaceae) in the central great barrier reef region.

BACKGROUND: Dinoflagellates are important primary producers, crucial in marine food webs. Toxic strains, however, are the main causative agents of non-bacterial seafood poisoning, a major concern for public health worldwide. Despite their importance, taxonomic uncertainty within many genera of dinoflagellates is still high. The genus Coolia includes potentially harmful species and the diversity within the genus is just starting to become apparent. METHODOLOGY/PRINCIPAL FINDINGS: In the current study, cultures were established from strains of Coolia spp. isolated from the central Great Barrier Reef (GBR). Cultures were identified based on thecal plate morphology and analyses of sequences (18S, ITS and 28S) from the nuclear rRNA operon. We report that the central GBR harbors a high diversity of Coolia species, including two species known to be capable of toxin production (C. tropicalis and C. malayensis), as well as the non-toxic C. canariensis. The strain of C. canariensis isolated from the GBR may in fact be a cryptic species, closely related but nevertheless phylogenetically distinct from the strain on which the holotype of C. canariensis was based. We also found evidence of the occurrence of a cryptic species morphologically very similar to both C. malayensis and C. monotis. The consequences of taxonomic confusion within the genus are discussed. CONCLUSION/SIGNIFICANCE: The central GBR region harbors a previously unreported high diversity of Coolia spp., including two species known to potentially produce toxins. The presence of a cryptic species of unknown toxicity highlights the importance of cryptic diversity within dinoflagellates.